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Abstract 

We study several phenomenologies of an Eg inspired extra U(l) model with 5*4 flavor symmetry. With the 
assignment of left-handed quarks and leptons to SVdoublet, SUSY flavor problem is softened. As the extra 
Higgs bosons are neutrinophilic, baryon number asymmetry in the universe is realized by leptogenesis without 
causing gravitino overproduction. We find that the allowed region for the lightest chargino mass is given by 
100-140 GeV, if the dark matter is a singlino dominated neutralino whose mass is about 36 GeV. 



* E-mail: yasu_daikoku@yahoo.co.jp 
^E-mail: hokada@kias.re.kr 



1 



1 Introduction 



Standard model (SM) is a successful theory of gauge interactions, however there are many unsolved puzzles in 
the Yukawa sectors. What do the Yukawa hierarchies of quarks and charged leptons mean? Why is the neutrino 
mass so small? Why does the generation exist? These questions give rise to the serious motivation to extend SM. 
Another important puzzle of SM is the existence of large hierarchy between electroweak scale Mw ~ 10 2 GeV 
and Planck scale Mp ~ 10 18 GcV. The elegant solution of this hierarchy problem is supersymmetry (SUSY)[1]. 
Recent discovery of the Higgs boson at the Large Hadron Collider (LHC) may suggest the existence of SUSY 
because the mass of Higgs boson; 125 — 126 GeV [2], is in good agreement with the SUSY prediction. Moreover, 
in the supersymmetric model, more information are provided for the Yukawa sectors. 

In the supersymmetric model, the Yukawa interactions are introduced in the form of supcrpotential. Therefore, 
to understand the structure of the Yukawa interaction, we have to understand the structure of superpotential. 
In the minimal supersymmetric standard model (MSSM), as the Higgs superfields H u and H D are vector-like 
under the SM gauge symmetry Gsm — SU(3) x SU(2) x U(l), we can introduce /x-term; 

HH U H D , (1) 

in supcrpotential. The natural size of parameter /x is O(Mp), however \i must be 0{Mw) to succeed in breaking 
electroweak gauge symmetry. This is so-called /x-problcm. The elegant solution of zx-problcm is to make Higgs 
superfields chiral under a new U(l)x gauge symmetry. Such a model is achieved based on ^-inspired extra U(l) 
model [3]. The new gauge symmetry replaces the //-term by trilinear term; 

XSH U H D , (2) 

which is converted into effective //-term when singlet S develops O(lTeV) vacuum expectation value (VEV) [4]. 
At the same time, the baryon and lcpton number violating terms in MSSM arc replaced by single G-intcractions; 

GQQ + G C U C D C + GU C E C + G C QL, (3) 

where G and G° are new colored superfields which must be introduced to cancel gauge anomaly. These terms 
induce very fast proton decay. To make proton stable, we must tune these trilinear coupling constants to be very 
small ~ O(10 -14 ), which gives rise to a new puzzle. 

The existence of small parameters in superpotential suggests that a new symmetry is hidden. As such a 
symmetry suppresses the Yukawa coupling of the first and the second generation of the quarks and the charged 
leptons, it should be flavor symmetry. We guess several properties that the flavor symmetry should have in 
order. At first, the flavor symmetry should be non-abelian and include triplet representations, which is the 
simple reason why three generations exist. At second, remembering that the quark and the charged-lepton 
masses are suppressed by SU(2)w gauge symmetry as the left-handed fermions are assigned to be doublet and 
the right-handed fermions are assigned to be singlet, the flavor symmetry should include doublets. In this case, 
if we assign the first and the second generation of the left-handed quarks and leptons to be doublets and the 
right-handed to be singlets, then suppression of Yukawa couplings is realized in the same manner as SU(2) W . 
At the same time, this assignment softens the SUSY-flavor problem because of the left-handed sfermion mass 
degeneracy. Finally, any products of the doublets should not include the triplets. In this case, we can forbid 
single G-intcractions when we assign G and G c to be triplets and the others to be doublets or singlets. As one 
of the candidates of the flavor symmetries which have the nature as above, we consider 54 [5]. In such a model, 
the generation structure is understood as a new system to stabilize proton [6] . 

In section 2, we introduce new symmetries and explain how to break them. In section 3, we discuss Higgs 
multiplcts. In section 4, we give order-of-magnitude estimates of the mass matrices of quarks and leptons and 
flavor changing processes. In section 5, we discuss cosmological aspects of our model. Finally, we give conclusions 
in section 6. 

2 Symmetry Breaking 

At first we introduce new symmetries and explain how to break these symmetries. The charge assignments of 
the superfields are also defined in this section. 
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2.1 Gauge symmetry 

We extend the gauge symmetry from Gsm to G32111 = Gsm X U(l)x X U(l)z, and add new superfields 
N c , S,G,G C which are embedded in 27 representation of Eq with quark, lepton superfields Q,U C , D c , L, E c 
and Higgs superfields H U ,H D . Where N c is right-handed neutrino (RHN), S is Gsm singlet and G,G C 
are colored Higgs. The two U(l)s are linear combinations of 11(1)^, U(l) x where E 6 D 5*0(10) x U{\)^ D 
SU(5) x f/(l)y x [7(1)^,, and their charges X and Z are given as follows 



X=^Q^ + \Q X , Z = -\q 4 , + ^Q x . (4) 

The charge assignments of the superfields are given in Table 1. To break U(l)z, we add new vector-like superfields 
$ c where < f )C is the same representation as RHN under the G32111 and its anti-representation $ is originated in 
27*. To discriminate between N c and $ c , we introduce Z-f symmetry and assign $ c , $ to be odd. The invariant 
superpotential under these symmetries is given by 



H / 32111 


= W + W S + W G + W*, 




yM 

h <S>$N C N C , 

M P 


(5) 




= Y U H U QU C + Y D H D QD C - 


h Y L H D LE C + Y N H U LN C - 


(6) 


W s 

w G 


= kSGG c + XSH u H D , 

= Y QQ GQQ + Y UD G C U C D C - 


V Y UE GU C E C + Y QL G C QL - 


_ y dn GD C N C , 


(7) 

(8) 




= M^ c + — Y*($$ c ) 2 , 
Alp 






(9) 



where unimportant higher dimensional terms are omitted. Since the interactions Ws drive squared mass of S to 
be negative through renormalization group equations (RGEs), spontaneous U(l)x symmetry breaking is realized 
and U(l)x gauge boson Z' acquires the mass 

m(Z') = 5V2g x (S) = 5%/2 ( Ap) ( S ) = °- 5255 ( S ) , ( 10 ) 

where the used value gx(Ms = ITeV) = 0.3641 is calculated based on the RGEs given in appendix A, and 
(H u - D ') <C (S) is assumed based on the experimental constraint 

m(Z') > 1.52TcV[7], (11) 

which imposes lower bound on VEV of S as 

(S) > 2892GcV. (12) 

To drive squared mass of $ c to be negative, we introduce 4th generation superfields and their anti- 

representations ,L4 and add new interaction 

W D Y LH ^ c Hl r L 4 . (13) 

To forbid the mixing between 4th generation and three generations, we introduce 4-th generation parity Z^ 
and assign all 4-th generation superfields to be odd. If M$ = in Wg>, then $, <!> c develop large VEVs along the 
D-flat direction of (<&) = ( ( E )C ) = V, U(l)z is broken and U(\)z gauge boson Z" acquires the mass 

m(Z") = 8g z V = 8 (l^9z\ V = 0.9202^, (14) 

where the used value gz(n = Mr) = 0.4365 is calculated by the same way as gx- We determine the values of 
two gauge couplings gx,9z by requiring three U(l) gauge coupling constants are unified at reduced Planck scale 
M P = 2.4 x 10 18 GeV as 

g Y {M P ) = gx(Mp) = g z {M P ). (15) 

In this paper we fix the VEV as 

V = Mi = 10 115 GeV. (16) 
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RHN obtains the mass 



M, 



M P 



10 4 " 5 GeV, 



(17) 



through the quartic term in Wq- 

After the gauge symmetry breaking, since the R-parity symmetry defined by 



R 



cxp 



20 



(3a: - 8y + 15 z) 



(18) 



remains unbroken, the lightest SUSY particle (LSP) is a promising candidate for cold dark matter. As we adopt 
the naming rule of superfields as the name of superfield is given by its R-parity even component, we call G, G c 
"colored Higgs". 

Before considering flavor symmetry, we should keep in mind following points. As the interaction Wg induces 
too fast proton decay, they must be strongly suppressed. The mass parameter M$ in W$ must be forbidden 
in order to break U(l)z symmetry. In W$, the contributions to flavor changing processes from the extra Higgs 
bosons must be suppressed [8]. 
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Table 1: G32111 assignment of superfields. Where the x, y and z are charges of U(l)x, U(l)y and U(l)z, and y 
is hypercharge. The charges of U{1)$ and U(l) x which are defined in Eq.(4) are also given. 



2.2 S4 flavor symmetry 

If we introduce S4 flavor symmetry and assign G,G C to be triplets, then Wq defined in Eq.(8) is forbidden. 
This is because any products of doublets and singlets of S4 do not contain triplets. The multiplication rules of 
representations of S4 are given in appendix B. Note that we assume full Eq symmetry does not realize at Planck 
scale, therefore there is no need to assign all superfields to the same flavor representations. In this model the 
generation number three is imprinted in G,G C . Therefore they may be called "G-Higgs" (generation number 
imprinted colored Higgs). 

Since the existence of G-Higgs which has life time longer than 0.1 second spoils the success of Big Ban 
nucleosynthesis (BBN)[9], S4 symmetry must be broken. Therefore we assign $ to be triplet and $ c to be 
doublet and singlet to forbid il/$$$ c . With this assignment, S4 symmetry is broken due to the VEV of <f> and 
the effective trilinear terms are induced by pentatic terms 

W NRG = i** c (GQQ + G C U C D C + GE C U C + G C LQ + GD C N C ) . (19) 

IVl p 

The size of effective coupling constants of these terms is given by 



m ($ c ) M 



M], " Mp " 10 . (20) 

This is the marginal size to satisfy the BBN constraint [10]. This relation gives the information about the RHN 
mass scale if the life time of G-Higgs is measured. 

The assignments of the other superfields are determined based on following criterion, (l)Thc quark and 
charged lepton mass matrices reproduce observed mass hierarchies and CKM and MNS matrices. (2)The third 
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generation Higgs H^,H^ are specified as MSSM Higgs and the first and second generation Higgs supcrficlds 
Hi 2 are neutrinophilic which are needed for successful lcptogenesis. 

To realize Yukawa hierarchies, we introduce gauge singlet and S4 doublet flavon superhcld Di and fix the 
VEV of A by 



V D = V\ ( D i) I 2 + I (D 2 ) | 2 = O.IMp = 2.4 x 10 17 GeV, (21) 
then the Yukawa coupling constants are expressed in the power of the parameter 

e = ^-=0.1, (22) 

which is realized by Zn symmetry. To drive the squared mass of flavon to be negative, we add 5-th and 6-th 
generation supcrficlds L^,q,D^ 6 as ^-doublets and their anti-representations L^ q^D^q and introduce trilinear 
terms as 

W 5 = Y^ID^DIDI + DIDD+D^DIDI-DIDI)] 

+ Y LL [D 1 {L 5 L 6 + L e L 5 ) + D 2 (L 5 L 5 -L 6 L 6 )], (23) 

where the mass scale of these fields is given by 

M Ls = Y DD V D = Y LL V D = eM P = 2.4 x 10 17 GeV. (24) 

We assign the 5th and 6-th generation supcrficlds to be Z\ -odd. The representation of all supcrficlds under the 
flavor symmetry is given in Table 2. The mass terms of 4-th generation fields are given by 

W A = Y LH ^ 3 L iH v + Y L {Dl +f* )2 L 4 L 4 + Y H {Dl ^ l)2 H¥ilY, (25) 

where 

M L4 = e 4 Y L M P = e 4 Y H M P = 2.2 x 10 14 GeV, (26) 
which realizes gauge coupling unification at Planck scale as 

gs(M P ) = g 2 (M P ). (27) 



2.3 SUSY breaking 

For the successful leptogenesis, the symmetry Z% x Z.^ must be broken softly. Therefore we assume these 
symmetries are broken in hidden sector and the effects are mediated to observable sectors by gravity. We 
introduce hidden sector superfields A, B + , Bi,-, B 2 _, C+, Ci,-, C 2 -, where their representations are given in 
Table 3. 

We construct O'Raifeartaigh model by these hidden sector superfields as follow [11] 
^hidden = ~ M ' '' A + m + B + C+ + mi-Bi-Ci- + m 2 -B 2 -0 2 - 



+ -\ + AC 2 + + -\ 1 -AC 2 _ + -\ 2 ^AC 2 _. (28) 
As the F-terms of hidden sector superfields given by 

Fa = -M 2 + h + C 2 + + h^Cf_ + ^Cl_ 7 (29) 

F B+ = m + C + , (30) 

F Bl _ = rm-Ci-, (31) 

F B2 _ = m 2 _C 2 -, (32) 

F c+ = m + B + + X + AC + , (33) 

F Cl _ = m-Bx- + Ai-ACi-, (34) 

F C2 _ = m-B 2 -+\ 2 _AC 2 ^, (35) 
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Table 2: S4 x Z 2 x x Z i7 x x Z% x ^2 assignment of superfields (Where the indices i and J of the 
^4 doublets runs % = 1, 2 and J = 5, 6 respectively, and the index a of the S4 triplets runs a = 1, 2, 3.) 
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Table 3: Z\ ] x x Z2 x C/(l)fl assignment of hidden sector superfields. All these superfields are trivial under 
the gauge symmetry G32111 and flavor symmetry S± x Z17 x Z-f" x Z^ x Z^ 5 ''. The observable sector superfields 
are Z^-evcn. 
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do not have the solution as 



(36) 

supersymmetry is spontaneously broken. The flavor symmetry x Z% is also broken. 

Since we assume the U(1)r symmetry is explicitly broken in the higher dimensional terms [12], soft SUSY 
breaking terms are induced by the interaction terms between observable sector and hidden sector as 



Csb = 



A 



w a w a 



M, 



■(cabcXaXbXc + ■ ■ •) 





Yh.c. 




. F 







A* A 

—^C a bX*Xb + h.c. 
1V1 1 ) 



(37) 



where the indices A, B, C runs the species of supcrficlds and the indices a, b, c runs generation numbers. Generally, 
as the coefficient matrices c a b are not unit matrices, large flavor changing processes are induced by the sfermion 
exchange. The explicit Z^ x Z 2 breaking terms are given by 

'B\B X - {D l X 1 +D 2 X 2 )*X 3 



C Z[ 2) B 



M 2 



M P 



h. 



D 



N R — 



= e m 2 BX (c(X 1 )*X 3 + S (X 2 )*X 3 ) + h.c. (X = H U ,H D ,S), 
^±^(c 2 jV| + c 3 NZ)*(N?) + h.c. 



Mi 



m 2 12 {N c 2 yN c 1+ mUN c 3 yNZ 



(38) 
(39) 



D 



2.4 S3 breaking 



The S 3 subgroup of 54 is broken by the VEV of ^-doublet flavon Di. Here we consider the direction of VEV. 
For the later convenience, we define the products of Di as follows, 



1 

T 



E 2 = D{ + D 2 , E 3 = 3D{D 2 - D% 
P 3 = D\ - 3D\D\, 

D 2 



V, = 



V 2 = 



2D 1 D 2 

12 n2 



Va 



-D 2 P 3 
DxP 3 



V s = 



-{D 2 ~D 2 )P 3 
2D X D 2 P 3 



(Di) = V D s = V D sixil 



D{-D 2 

and the VEVs of each components of Di as 

(£>i) = V D c = V D cos 

Generally, the superpotential of Di is written in the form of polynomial in E 2 ,E 3 , P 3 as 

M P 4 W D = a x E\E 3 + a 2 E\E\ + a 3 E*P 2 E 3 + a A E 2 El + a 5 E 2 P 2 El + a 6 E 2 P^E 3 
Substituting the VEVs given in Eq. (43) to the flavon potential, we get 



(40) 
(41) 

(42) 



(43) 
(44) 



V(Vb, 9) = | -A[a' lS3 + a 2 s 3 3 + a' 3 c 2 3 s 3 + a' 4 s 5 3 + a' 5 4s 3 3 + a' 6 c 4 3 s 3 ]Vl, 
where 



(Ye.) 
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V F + m 2 V 2 , (45) 



S3 = sin 39, c 3 = cos 39, 
and Vf is F-tcrm contribution. As this potential is polynomial in S3, the stationary condition 
dV(V D ,9 



09 



\ II 1 11 1 II 2 1 11 3 11 4 1 // 5 1 11 7 , 11 91 n 

c 3 [a + a 1 s 3 + a 2 s 3 + a 3 s 3 + a 4 s 3 + a 5 s 3 + a 7 s 3 + a 9 s 3 \ = 0, 



gives parameter independent solution 

c 3 = 0, 

and parameter dependent solution 

a'o + a"s 3 + a 2 s 2 3 + a 3 sl + a^s 3 + etgSg + a"s\ + a'^sl = 0. 



(46) 

(47) 

(48) 
(49) 
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Which solution of two is selected for the global minimum is depends on the parameters in potential. Since the 
solution Eq.(48) gives wrong prediction such as masslcss up-quark and electron, we assume the solution Eq.(49) 
corresponds to the global minimum. In this paper we assume (Di) are real without any reason, which is important 
in considering CP violation in section 4. 

The scale of Vjj is determined by the minimum condition 



1 dV(V D ,9) _ m2 Vj3Q 
V D dV D M P 



D 



2s 



0, (50) 



as 

Vd_ [H\ 1/15 ( 10 3 GcV \ 1/15 
M P ~ \M P ) ~ I 10 18 GcV/ 



10" 1 , (51) 



which agrees with Eq.(22). In this paper we sometimes write SUSY breaking scalar squared mass parameters as 
m 2 for simplicity and assume m ~ O(TeV). 

2.5 S4 breaking 

The superpotential of gauge non-singlet flavons $, $ c is given by 

W * = £( $ 3) 2 [*l+*2 + *3] + ^[(^) 2 + ( < I>2) 2 ][ < I'l+ < i>2 + < i ) 3] 

+ 2L[2Vm»?«I>§(#1 - $ 2 ) + (M) 2 - (<^) 2 )($ 2 + $ 2 - 2$ 2 )] 

+ ^#§[\/3$;($i-$i) + #§($l + $i-2*?)]. (52) 

Since the first term in Eq.(25) drives the squared mass of $3 to be negative through RGEs, these flavons develop 
VEVs along the D-flat direction as follows 

m) = = 0, ($0 = ($ 2 ) = <$ 3 ) = M = ii, (53) 

where ^-symmetry is unbroken in this vacuum. The scale of V is determined by the minimum condition 

1W( * ) -- a "y|»^ = J (54) 



V <9$ Mp 



as 



V_ \m\ 10 3 GcV 7 

A/ P ~ Y 1^*1-^ ~ V (0-l)10 18 GeV ~ ' 1 ' 

which agrees with Eq.(16). In this paper we define the size of O(l) coefficient as 0.1 <Y X < 1.0. 
Note that there are S3 breaking corrections in the potential of $, <& c as follows 



2 ( (pm + p 2 mY(Din) + (gig! - ggjircggj^ 

|2£ 2 $i| 2 + |(V3£>i + £ 2 )$ 2 | 2 + \(VSD 1 - £ 2 )$ 3 | 2 ' 



2 \ w 



e 2 m 2 [s 2 (^)*$ c 3 + c 2 ($ c 2 )*$ c 3 + h.c.] 

e 2 m 2 [4s 2 |$ 1 | 2 + {Vic + s) 2 |$ 2 | 2 + (Vic - s) 2 |$ 3 | 2 ] + ■ ■ ■ , (56) 



the direction given in Eq.(53) is modified as follows 

m~m)~0(e 2 )V, m) = (l + 0(e 2 ))V, ($ Q ) = (1 + 0( e 2 ))^=. 
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3 Higgs Sector 



Based on the set up given in section 2, we discuss about phenomenology of our model. In this section, we consider 
Higgs doublet multiplets H^,H^ and singlet multiplets S a - 

3.1 Higgs sector 

The superpotential of Higgs sector is given by 



Ml Si 



A" 

^S 3 [{2D X D 2 ){HYH? + H%H?) + {D\ - Dl)(H?H? H% H?)] 
mi 

S^D^ + D 2 H\ I )(D 1 Hi ) + D 2 H?) + (Drf? - D 2 H^){D l H^ ) - D 2 H?)\ 



+ \ 3 S 3 H 3 H 3 + X^H 3 (SiH 1 + S 2 H 2 ) + \5(SiH 1 + S 2 H 2 )H 3 

+ kS 3 (G 1 Gl + G 2 G^ + G 3 Gl). (58) 
For simplicity we assume 

X[ = X" = X"' = X'l" = 0, (Ax) e 2 = A!. (59) 

The coupling k and A3 drive the squared mass of S3 to be negative. 
Omitting 0(e)-terms, Higgs potential is given by 



V = m^(|^| 2 + |^| 2 ) + m^|^| 2 + m^4|Ff| 2 + |i/f| 2 )+m^ f |^ 
+ mldStf + lS^+mljSsf 

{X 3 A 3 S 3 H^H 3 D + XiMHYiSxH? + S 2 /ff ) + )*M{SiH? + S 2 H^)Hi> + h.c.} 
+ J X 3 H 3 H 3 I + \X&H 3 + X^H^ Hi \ + \\±H 3 H 2 + X§H 3 H 2 \ 
+ \X 3 S 3 H° + A 4 (5iFf + S 2 i/ 2 D )| 2 + IXsHgSA 2 + lAsi^l 2 
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+ \X3S 3 H 3 + Xr,(SiH 1 +S 2 H 2 )\ +\\4,H 3 Si\ +\X4H 3 S 2 \ 
+ £ )Wtf + (H??a A H?] 2 + \g 2 Y [\H^\ 2 - |i^| 2 ] 2 

A=l 

+ \gl[x u \HV\ 2 + x d \H°\ 2 +x s \S a \ 2 } 2 + V lAoo ^ (60) 

where Vj_i 00 p is 1-loop corrections from Q 3 ,U 3 ,G a ,G c a . The VEVs of H 3 ,H 3 ,S 3 trigger off gauge symmetry 

(2) 

breaking at low energy scale. Z\ -breaking terms 

V FB = em 2 BU (H^c + H 2 s)*H 3 + em 2 BU (H°c + H 2 s)* H 3 + em 2 BS (S lC + S 2 s)*S 3 + h.c, (61) 

enforce SVdoublets developing VEVs as follows 

2 \ 

BX \ i v \ 1 v \ „ / c " l BX 



(Xx)~c )(X 3 ), (X 2 )~s — p )(X 3 ), X = H U ,H U ,S. (62) 

V m x J V m x J 

Due to the symmetry, the Yukawa couplings between H 3 and N c are forbidden and neutrino Dirac mass is 
not induced. To give neutrino Dirac mass, we assume the size of VEV of Hf is given by 

(if^-e^^-lGeV, (63) 
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(2) 

and put the Z\ breaking parameters as follows 

m 2 BU ~ m% D ~ m% s ~ «™|uSY' ( 64 ) 

by hand. The suppression factor O(e) may be induced by the running based on RGEs, because off diagonal 
elements of scalar squared mass matrix do not receive the contributions from gaugino mass parameters which 
tend to make scalar squared mass larger at low energy scale. 
We use the notation of VEVs as follows 

(Hf) = {c,s)v u , (H^) = v' u , (H t D ) = (c,s)v d , (H 3 D )=v' d , = (c,s)v s , (S 3 )=v' s , (65) 

where we fix the values by 

v' u = 150.7, v' d = 87.0, v' s = 4000, v E w = J(K) 2 + K) 2 = 174 ( GcV )> tan /? = ^ = tan ~. (66) 

In this paper, we neglect the contributions from v u ^, s except for neutrino sector. With this approximation, the 
potential minimum conditions are given as follows, 

= v u my = ^ " XM {v ' d/v ' u) + xl{v ' d? + A ' K)2 + + gl)[Kf {v ' d?] 

2.^[-2K) 2 -3W) 2 + 5K) 2 ] + ^^-^P, (67) 

= V d M^ = m k ^MK/v' d ) + \IK) 2 + x 2 M) 2 \{g 2 y + gl)\K) 2 - K) 2 ] 

3 5 2 [-2«) 2 - 3(v' d ) 2 + 5(v>) 2 } + (68) 

o = ^r|^ = m ^-\MM<) + xl{v l J 2 + \l{v' d ) 2 + ^ (69) 

v s ob 3 

where the 1-loop contribution is neglected in Eq.(69), which is unimportant. These equations give the boundary 
conditions for m 2 HU ,m 2 HD ,m S3 at SUSY breaking scale Ms — 10 3 GcV in solving RGEs. 
The mass matrices of heavy Higgs bosons are given as follows 

/ \ 3 A 3 v' s v' d K -X 3 A 3 v' s \ 

M|(CPeven) ~ -X 3 A 3 v' s X 3 A 3 v' s v'Jv' d (70) 

V 50. 9 2 K) 2 J 

M 2 (CP even) = M 2 (CP odd) ~ diag (m 2 HU - 10g 2 x {v' s )\™ 2 HD - 15g 2 M) 2 ,™ 2 + 25 5 2 K) 2 ) (71) 

/ v'J< 1 \ 

M|(CP-odd) ~ X 3 A 3 v' s 1 v'Jv' d (72) 

\ / 

Af 2 (charged) ~ X 3 A 3 v' s ( V ' d/ < ^ ) (73) 

M 2 (charged) ~ diag (m 2 H u - 10g 2 x (v' s ) 2 , m 2 H o - 15g 2 x (v' s f) , (74) 

where only O(lTeV) terms are considered and notation H U H D = (H U )°(H D )° - (H U )+(H D )- is used. In this 
approximation, generation mixing terms are negligible. The third generation mass matrices are diagonalized as 
follows 

Mf(CPeven) = diag (0, ^r^, 50 ff 2 (^) 2 ) , (75) 
Mf(CPodd) = diag fo.^M.o), (76) 

Mf (charged) = diag (o, ^^f) , (77) 

where the zero eigenvalue in CP even Higgs bosons corresponds to lightest neutral CP even Higgs boson and the 
other zero eigenvalues are Nambu-Goldstonc modes absorbed into gauge bosons. 



10 



3.2 Lightest neutral CP even Higgs boson 

To calculate the mass of the lightest neutral CP even Higgs boson, we diagonalize 2x2 sub-matrix given by 



A/3 (even) 



hid, 3 



l uu,3 
l ud,3 



'ud, 3 
1 dd, 3 



m uu,3 — 



L ud,3 



-A 3 A 3 < + 2Xlv' u v' d - l -(gl + gl)v' u v' d + I2c? x v' u v d + L l ~ lo01? 



2 dv'dv' 



hA3(v' s v'Jv' d ) + l(g 2 Y + g 2 2 )(v' d f 



ld 2 ^ 



l-loop 1 ^^i-loop 



d(v'} 2 



2v' dv' 



(78) 
(79) 
(80) 
(81) 



J di " u d ""d 

where 0{vew) terms are included. We evaluate the l-loop contributions from top, stop, G-Higgs and G-higgsino 



^l-loop ~ 7^2 Str 



AP In 



Mr 
A 2 



(82) 



at the renormalization point A = Ms [13] [14], where the mass eigenvalues are given by 



m T± =M T + (Y 3 u v' u y±R T , 
Ml = Hm^+my+SggK ) 2 , 
Rt = ^(AA/ 2 ) 2 + (Y 3 U X T )\ 
AA/ 2 = 



m G ± = M G ± Rg, 



U m Q 3 



l u 3 ) 



X T = \3v' s v' d - A"v' u , 
m t = Y 3 u v' u , 



Ml 

R G = V(AA/ 2 ) 2 + {kX G f, 
AM G = i (m 2 G - rap) + §3 2 (ws) 2 , 
X G = hi'Li'd - A k v' s , 
m„ = kv' 



(83) 



where we neglect 0(vew) terms in D-term contributions. The mass matrices of stop and G-Higgs are given in 
following sections (see Eq.(171) and Eq.(236)). By the rotation with 



V, 



11 



VEW 



the O(lTeV) term is eliminated from off-diagonal element of Eq.(78) and we get 



= (V£M 2 V H )n = (m 2 )o + (m 2 h ) T + (m 2 )c 



(A3 sin pf 



9y + 9-2 



cos 2 2/3 + 2 5 2 (2 sin 2 /? + 3 cos 2 (if 



J EW 



(84) 

(85) 
(86) 



(ml) 



U\2 



T — 



3(1?) 
16ir 2 v 2 EW 



(Y?X T ) 2 
R\ 



+ 2(1?W J )> 



(Y 3 u v u y 



In- 



2(Y 3 v'J X T _ _ 3 + {Y u v ^ <Y a X T )' 



( m h)c = 



Rt 

9fc 2 (A 3 <^) 2 
Sk 2 v 2 ew 



2R T 



In- 



(87) 



-2 



(AA/ 2 ) 2 



r g 



R% 



In 



In- 



G 



A 2 



In. 



A 2 



If we fix the parameters at Ms as given in Table 6, we get 

82.7, v /(m 2 ) T -94.2, y / K)^ = 9.2 (GcV). 



m ft = 125.7, 



(89) 



The l-loop contribution is dominated by stop and top contributions, this is because we put k small (k = 0.5) 
to intend the mass values of the particles in the loops are within the testable range of LHC at y/s = 14TeV as 
follows 



m T+ = 1882, m T _ = 1178, m G+ = 3908, m G _ = 1737, 



2000 (GeV). 



(90) 



The value of A3 = 0.37 is tuned to realize observed Higgs mass which is mainly controlled by this parameter 
through (\3VEW sin/3) 2 and Xt for fixed v' s and A 3 (= A t ). 
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3.3 Chargino and neutralino 

At next we consider the higgsinos and the singlinos. The mass matrix of the charged higgsinos is given by 



C 



c 



((hV) + ,(h%) + ,(hV)+) 



Xiv' s X 5 v s c 
Xiv' s X 5 v s s 
X^VgC X4V s s X 3 v' s 




+ h.c. 



(91) 



Since the (3,3) element is much larger than the other 0(e 2 ) elements, the first and second generation higgsinos 
decouples and have the same mass X\v' s . With the gaugino interaction given as follows 



C 



gaugino 



1 — X A X A 



92 *>2 T 2 + ^9yX y - 2g x X x 



A=l 
3 

92 ^ 

L A=l 

1 



3 1 

^ ^ 2 T 2 ~ 7^9yXy - 3g x Xx 



h? -iV2(S a )^[5g x X x ] Sa 



2 -M 2 A^ \* - -M Y X Y X Y - -M x XxXx + h.c, 
the third generation charged higgsino mixes with wino and the mass matrix is given by 



w 



,w 
—i 

71 



A 3< g 2 v' u 



92V d 



+ 2 



(h 



D\- 
3 I 

W~ 



+ h.c. 



(Xl T iX 



(92) 

(93) 
(94) 



The mass eigenvalues of charginos are given by 

1 



M 



xf 



[M s f+gl{v' u Y+gl{v' d y + M 



± \/\[(hv' s ) 2 + ( 5 2<) 2 - (92v' d f - M|]2 + (X 3 g 2 v' s v' d + M 2 g 2 v' u f 
M(xf) = Xiv'., 



(95) 
(96) 



where X3 is almost third generation higgsino and Xw ^ s almost wino. 

The mass matrix of the ncutralinos is divided into two 3x3 matrices and one 6x6 matrix as follows 







Xiv' X 5 v d 



( 



M 



T 

Xo 



i=i,2 \ X 5 v' d A 4 < 

X s v' s X 3 v' d g Y v'jV2 

X 3 v' s X 3 v' u ~9Yv' d /V2 

X 3 v' d A 3 < 

g Y v' u /V2 -g Y v' d /V2 -My 

-g 2 v' u /V2 g 2 v' d /V2 

V -2g x v' u V2 -3g x v' d V2 5g x v' s V2 

((h^) ,(hi) ,s 3 ,iX Y ,iXliX x ). 




- 7;Xo M xXo, 



-g 2 v' u /V2 ~2g x v' u V2 \ 

g 2 v' d /V2 -3g x v' d V2 

5g x v' s V2 



-M 2 

-M x J 



(97) 



(98) 



(99) 



The mass eigenvalues of these mass matrices are given in Table 7. The common mass of two LSPs is given by 
the smallest eigenvalue of 3 x 3 matrix given in Eq.(97). Note that the LEP bound for chargino 



Xiv' s > lOOGeV [15], 



(100) 



must be satisfied. Requiring the coupling constants X^§ do not blow up in fi < Mp, we put upper bound for 
them as A4 < 0.5, A5 < 0.7, then the rough estimation of LSP mass is given by 



M(x o } „ 2(XW u )(X 5 v' d ) < 9000GeV 2 



(101) 
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where Xi 1 i s almost singlino. To realize density parameter of dark matter QcDMh = 0.11, we must tune 
M(Xi i) ~ 30 — 35GeV to enhance annihilation cross section. This condition gives upper bound as 

\!v' s < 300GeV. (102) 

This constraint is not consistent with the lower bound as follows 

m v ± > 295GeV[16], m v ± > 330GeV[17l. (103) 

Xi Xi 

Therefore we assume the lightest chargino mass is in the region 

100 < X lV ' s < 140 (GeV), (104) 

in which 3-lepton emission is suppressed due to the small mass difference between chargino and neutralino 
compared with mz- 

Note that bino-like neutralino can decay into Higgs boson and LSP through the 0(e 2 ) mixing of Higgs bosons 
by the interaction 

C D -iV2(H?)° (\9yXy) (h?)° ~ 0(e 2 )(HV)°\ Y (h?)°. (105) 

4 Quark and Lepton Sector 

In this section, we consider the quark and lepton sector and test our model by observed values given as follows, 
running masses of quarks and charged leptons at fi — Ms = ITeV [18] 

m u = 1.10t^(MeV), m c = 532 ± 74(McV), m t = 150.7 ± 3.4(GcV), 

m d = 2.50li;^(MeV), m s = 47±i3(MeV), m b = 2.43 ± 0.08(GeV), (106) 
m e = 0.4959(MeV), m M = 104.7(MeV), m T = 1780(MeV), 

CKM matrix elements [19] 

\V ud \ = 0.97427, |K S | - 0.22534, \V ub \ = 0.00351, 

\V cd \ = 0.22520, \V CS \ = 0.97344, jv^l = 0.0412, (107) 

\Vtd\ = 0.00867, |VJ,| = 0.0404, \V tb \ = 0.999146, 

and neutrino masses and MNS mixing angles [19] 

Ami, = < " < = (7-^-°of 6 ) x 10" 5 (eV 2 ), (108) 
Am 2 2 = |m 2 3 -m 2 2 | = (2.35^ 2 9 )xl0" 3 (eV 2 ), (109) 




C12C13 S12C13 

VmNS = -S12C23 - Ci 2 S23Sl3e l5 C12C23 - Si 2 S23Sl3e l5 S23C13 

S12S23 - Ci 2 C23Sl3e*' 5 -C12S23 - Si 2 C23Sl 3 e 1 ' 5 C23C13 

sin 2 9 12 = 0.306±g;81|, sin 2 8 23 = 0.42±g;8|, sin 2 8 13 = 0.02lt°;°°£. ( n 0) 
After that we estimate the flavor changing process induced by sfermion exchange. 

4.1 Quark sector 

The superpotential of up-type quark sector is given by 

Wu = Y 3 U H^Q 3 U^ + e 2 Y 2 u HV [Q lS2 + Q 2 c 2 ]U c 3 + e 3 Y^H 3 [Q lC + Q 2 s]U^ + e*Y?HgQ 3 U? 

+ e 6 {Y^H^[Q lS2 + Q 2 c 2 \U{ + Y 6 u H^[s 3 (Q lC + Q 2 s)]U^ - Y 7 U H^[c 3 (Q lS - Q 2 c)}U<} , (111) 

from which we get up-type quark mass matrix as 

e 6 (Y 5 u s 2 + Y 6 u s 3 c ~ Y 7 u c 3 s) e 3 Y 4 u c e 2 Y 2 u s 2 ' 
M„=| e e (Y 5 u c 2 +Y 6 u s 3S + Y 7 u c 3 c) e 3 Y^s e 2 Y 2 u c 2 \ v' u = \ e 6 e 3 e 2 \ v' u . (112) 
e 4 ^ Y 3 U J \ e 4 1 



13 



Note that there is dangerous VEV direction such as 8 = ^. In this direction the matrix given in Eq.(112) is given 
by 

/ e 6 (YU + Y 6 u )c e 3 Y^c e 2 Y 2 u c \ 
M u =\ e e (Y 5 u + Y 6 u )s e 3 Y" s e 2 Y 2 u s \v' u , (113) 
V Y? J 

which has zero eigenvalue. In the same way, the down type quark masses are given by 

W D = e 4 {Y 2 D Hi ) [s 3 (Q 1 c + Q 2 s)]D c 3 + Yi > H 3 u {Q 1 s 2 + Q 2 c 2 )D c 3 + Y g D H 3 [c 3 (-Q 1 s + Q 2 c)]D c 3 
+ Y s D H^[a 3 (-Qia + Q 2 c)]D c 2 + YfHfi-Qia + Q 2 s 2 )D c 2 + Y^H^c^Qxc + Q 2 s)]D c 2 ) 
+ e 5 {YfHglQxc + Q 2 s]Dl + Y 7 D H^s^Q^ + Q 2 c 2 )\D\ 

+ Y£H 3 D [c 3 (-Q 1 c 2 + Q2S 2 )m} + e% D H»Q 3 D c 3 + e 3 Y 1 D S3 H?Q 3 D c 1 , (114) 
from which we get down-type quark mass matrix as follows 

Y°c 3 c 2 ) e 4 (-Y 5 D s 3 s-Y 6 D c 2 +Y°c 3 c) e\Y 2 D s 3 c + Y A D s 2 - Y 9 D c 3 s) 
■Y£c 3 s 2 ) e 4 (Y 5 D s 3 c + Y 6 D S2 + Y»c 3S ) e 4 (Y 2 D s 3 s + Y 4 D c 2 + Y 9 D c 3 c) 

e 2 Y 3 D 

(115) 



M d /v' d 






( 


Y 7 D s 3 s 2 


+ Y 8 D c 




^ 7 D s 3 c 2 


+ Y 8 D s 
e 4 Y, D 










e 4 e 4 










(2 


e 4 e 4 


) 




e 2 





The effects of flavor violation appear not only in supcrpotcntial but also in Kahlcr potential as follows 

k{u c ) = \utf + \uz\ 2 + \ut\ 2 - 



K{D C ) = \D\\ 2 + \D C 2 \ 2 + |D C ' 2 




(C/ 3 T)| £ 3 1 e 5 



M 3 M% 



1 


e 3 




e 3 


1 


?)(!) 


e 3 







((Dtr, (!?§)*, (Z>f)*) e 3 1 e 6 £>2 , (117) 



K{Q) = {\QiY + \Q 2 \ 2 ) + \Q 3 \ 2 + \ v Z' ^° + ' \ r2 ' +--- + h.c 




(Ql,Ql,Ql) | ^ i e 2 ) I Q2 | , (118) 

where dot in X ■ Y means inner product of two S^-doublets X, Y. Therefore, a superfield redefinition has to be 
performed in order to get canonical kinetic terms as follows [20] 

/ {my \ / i e 3 e 4 \ 

= V K (U) my , V K {U) = e 3 1 e 5 , (119) 

V (my J \e 4 e 5 1 J 

( (my \ / i e 3 e 3 \ 

= V K (D) {D c 2 y , V K (P) = e 3 1 e 6 , (120) 

V m' J \ e 3 e e ! y 

= VOf(Q) L MQ) = e 2 1 e 2 , (121) 
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by which the mass matrices given above are transformed into 



M' u = V£{Q) \ e 6 e 3 e 2 \v' u V K {U)=\ e 6 e 3 e 2 \v' u , (122) 
4 1 / \ e 4 e 5 1 








4 4 \ / 5 4 4 

rT(rt\ I ^5 ^4 ^4 I Jt/ /r>\ _ I J> A ^4 I „,/ 



M' d = V£(Q) 6 5 e 4 e 4 ^ (0) = e 5 e 4 e 4 t£. (123) 



,-2 



These matrices are diagonalizcd as follows 



1 e 3 e 

M' u = L u M^Rl = ( 1 1 e 2 ] ( e 3 1 < ( e 3 1 e 5 I , (124) 



1 


1 










1 


1 


5 ) ( : 


c 3 





r 2 


c 2 







1 


1 


1 










1 


1 


? ) ( : 


e 4 





c 2 


r 2 







e 2 



,4 £ 5 



= L dM f^R\=[ i i e 2 I [ e 4 U( e 3 1 e 4 ). (125) 



e e 4 1 



Therefore Yukawa hierarchies arc given by 



F U (M F ) = e 6 , F C (M P ) = e 3 , F t (M P ) = Y^{M P ) = 1, 

r d (M P )-e 5 , F s (M P )=e 4 , n(A/ P ) = e 2 . (126) 



On the other hand, observed values Eq.(106) give 



y,(m p) - ±(±l^)= 1Al °, (127) 



5.1 V 150.7 

F t (Mp) = 0.28, (129) 
M M P) . (130, 



7.2 V 87 



7.2 V 87 

F d (Mp) = -L = 0.39e 2 , (132) 



7.2 V 87 . 

which give good agreement with Eq.(126). Where the renormalization factors 



./^=5.1, ,/^ = T.2, (133) 
y a u (Mp) y a d (Mp) 

and F(Afp) = 0.28 arc calculated based on RGEs given in appendix A. CKM matrix is given by 

/ 1 1 e 

Vckm = L\L d = 1 1 e 2 | , (134) 
V e 2 e 2 1 

;h Eq.(107). 

(2) 

Note that the breaking induces generation mixing in Higgs bosons then Yukawa interactions arc modified 



which is consistent with Eq.(107) 

Note t 
as follows 



-C = Y% {H$ + t 2 H l ( + e 2 H^) qi u c 3 + Y° (H° + e 2 H» + e 2 H 2 D ) qi d c y (135) 

Since these Yukawa coupling matrices are diagonalized in the basis that the quark mass matrices are diagonalized, 
the extra Higgs boson exchange do not contribute to the flavor changing processes. 
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4.2 Lepton sector 

With the straightforward calculation, the mass matrices of lepton sector are given as follows. From the superpo- 
tentials 



We 



W n 



e 3 H 2 (la i L 2 , L 3 



e 5 (Y E c 

'E 



Y 8 E s 3 s 2 - Y E c 3 c 2 ) 



e 5 (Y E s + Y E s 3 c 2 + Y E c 3 s 2 ) 
e 5 Y E s 3 

Yfss + Y 4 N cs 2 + ■■■ 
-Y 2 N c 3 + Y 4 N cc 2 + ■■■ 



e 3 Y E c 
e 3 Y 5 E s 



e 2 Y 4 E s 2 
e 2 Y E c 2 



e 2 Y E 



Y 3 N c 



e 3 Y E s 3 

Y 5 N s 3 + Y 8 N cs 2 
-Y 6 N c 3 + Y 8 N cc 2 
Y 7 N c 



Y 4 N ss 2 



vN„ 
Y 2 c 3 

Yfss + Y 4 N sc 2 



Y 6 N c 3 + Y s N ss 2 

rN , 



Y*s 



Y»s 3 + Y»sc 2 



Y 7 N s 




($2 + $2 + $2) [YgNZN? + Y 2 N 2 N^N* + Y 3 N 3 N'Nc + Y 2 N 3 N C 2 N : 



we get original mass matrices as follows 



Mr 




M 



M 




Redefining the Kahlcr potential given by 



K(E C 



K(L) 




f {E 2 El)*El 


(E 3 Ef)*EI 


\ Ml 


Ml 


1 e 2 £ 3 \ 




£ 2 1 e 3 


(§)• 


e 3 e 3 I J 





(y 2 E c 2 y ■ (ViE%) 



Ml 



h.c. 



L 2 D 1 y(D 1 L 3 ) + {L 1 D 1 -L 2 D 2 y{D 2 L 3 ) , \L-V X \ 



K{N C ) = 

= ((iVf)*,(iVf)*,(iV 3 c )*) 
by the superfields redefinition as 

= V K {E) 
= V K {L) 




M 2 P 



h.c. 



h.c. 




= V K (N) 




V k {E) = 



V K (L) = 



V K (N) = 



1 


e 2 


c' 


c 2 


1 


e 3 


e 3 


e 3 


1 








e 2 


e 2 




1 


r 2 


) 


c 2 


1 




1 











1 


r 2 





e 2 


1 



the modified mass matrices are given by 
M> = Vg(L) 



v' d V K {E) 



(136) 



(137) 
(138) 

(139) 



(140) 



(141) 



(142) 

(143) 
(144) 
(145) 

(146) 
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1 1 \ / 1 1 \ 

M' D = V£(L) | 1 1 e 3 v u V K (N) = 1 1 e 3 w„, (147) 

11/ \ 1 1 / 

1 \ 2 / 1 \ 2 

M> M = V-(N) | 1 1 J —V K (N) = ^ 1 1 j — . (148) 

The mixing matrices of charged leptons are given by 

M' e = L e M^R\ = ( 1 1 1 e 3 j v' d ( e 2 1 e 

e 2 / \ e 3 e 1 

The Yukawa hierarchy of charged leptons gives good agreement with the experimental values given by 

y AMp) _ 1(2^2^^0.30^ (150) 



1.9 V 87 
J_ / 1.78 



^r(Afp) = J- ( ^f) = 1.08e 2 , (152) 



where the used value of rcnormalization factor 



'ae(Ms) 

ffle(Mp) 



1.9 (153) 



is calculated based on RGEs given in appendix A. 
The neutrino seesaw mass matrix is given by 



1 1 1 
1 1 1 
1 1 1 



m v = = O(0.01cV), Mr = (154) 



which has one zero eigenvalue because one RHN n\ does not couple to left-handed leptons. Therefore mixing 
matrix and mass eigenvalues are given as follows 



L V M V L V = diag(m ! y 1 , m V2 , m„ 3 ) , (155) 
1 1 1 

L v = 111 (156) 
1 1 1 



m U2 = sjrri^ = 0.87 x 10~ 2 , m V3 ^ v/mj^ = 4.8 x 10~ 2 (cV), (157) 



and MNS matrix is given by 



1 1 1 

Vmns = L\L V = I 1 1 1 I . (158) 



1 1 1 

With the recent experimental value of | sin #13! ~ 0.14 [21], MNS matrix is given by 



0.64 0.55 0.14 

V MN s = \ 0.42 0.64 0.65 ] , (159) 
0.36 0.55 0.76 



which is consistent with Eq.(158). 
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4.3 Squark and slepton sector 

Sfcrmion mass matrices are given as follows 

-C 



D 

M 2 {Q) = 

M 2 (U C ) = 

M 2 {D C ) = 

M 2 {L) = 

M 2 {E C ) = 

A(U) = 

A(D) = 

A(E) = 
V F 
V d 



X* a [M 2 (X)] ab X b 



mQ a A{U) ab U c b 

„,2 



H^Q a A{D) ab D c b + H»L a A{E) ab E c b 



m 2 e 2 



2 2 

m z e z 



m 2 e 2 



2 2 

m z e 



m 2 e 2 



m 2 e 2 



m z e J 



2 2 

m z e 



2 2 

me 

2 2 

m z e z 



m 2 e 3 
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e 2 Y u A u 
e 2 Y u A u 

YfA" 



E l 

e 3 Y u A u 
e 3 Y u A u 



e 5yD^D e 4 Y D A D e 4 Y D A D 

e 5 Y D A D e A Y D A D e 4 Y D A D 


e 3 Y E A E 

e 3 Y E A E 
e 3 Y E A E 



e e Y u A u 

e e Y u A u 
e i Y u A u 



e 3 Y D A D 



e 5 Y E A 



e 5 Y E A 



e 5 Y E A 



2 Y D A D 



e 2 Y E A 
e 2 Y E A 
e 2 Y E A 



= \Y 3 U HVQ 3 



\Y 3 U H^U^\ 2 



\y 3 u Q^ 



\ 3 S3H 3 \ 2 , 



1 2 

2 9 * 



5|5 3 | 2 



^(|Q a | 2 + |C/ c | 2 + 2|^| 2 + 2|L a | 2 + |^ 



h.c.] + V F + V D , (160) 
(161) 

(162) 

(163) 

(164) 

(165) 

(166) 

(167) 

(168) 
(169) 
(170) 



a=l 



where m = O(10 3 GeV) and the contributions from F-terms are neglected except for top- Yukawa contributions 
and the contributions from D-terms arc neglected except for the contributions from S3. 

After the redefinition of Kahlcr potential and the diagonalization of Yukawa matrices, sfcrmion masses arc 
given as follows 



I0g 2 x (v' s ) 2 ) K| 2 + EW + ^K) 2 )IQ.f 



i=l 



i = l 



K 3 + ^ 2 M) 2 ) \d 3 \ 2 + J2 K + wglWs) 2 ) m 



i=l 
3 



« + I0g 2 x (v' s ) 2 ) \L 3 \ 2 + £ ( m k + ^ 2 M) 2 ) \K\ 2 



a=l 

m 2 Qs + (Y 3 u v' u ) 2 + 5gl(v' s ) 2 Y 3 u \ 3 v>v> d - A u 3 Y 3 u v' u 



U 3 

my 

m 2 U* a {5 u LL ) ab U b + m 2 D* a (5^ L ) ab D b + m 2 {U c )* a {5 u RR ) ab U^ + m 2 (D c ):(Sji R ) ab D c b 



Y 3 u X 3 v' s v' d ~ A u 3 Y 3 u v' u m 2 Ua + {Y 3 u v> u ) 2 + 5g 2 x {v> s ) 2 
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°LL 



m 2 E*M L ) ab E b + m 2 N* a {8» L ) ab N b + m 2 (E c y a (5% R ) ab E c b 

m 2 {U a {8 u LR ) ab U c b + D a (8l R ) ab D c b + E a (Sf R ) ab E c b + h.c.} , (171) 

1 / e 2 e 2 \ 

^[LiV^Q)M 2 (Q)V K (Q)L u ] oSdi&Kon&1 = ( e 2 e 2 I , (172) 



<*£z = i[44(Q)M 2 (Q)^(Q)id off diagona i = I e 2 e 2 J , (173) 



e 2 e 2 
e 2 




1 / 1 1 \ 

5h = -s [44 (L)M 2 (L)V K (L)L e ] oS di j = 10 1, (174) 

m \ 1 1 J 

1 / £ 2 e 4 \ 

— [RM(U)M 2 (U)V K (U)R u ] oSdisLROnal ^ \ e 2 e 5 , (175) 



°RR 



°RR 



e 5 

e 3 e 

1 "?t v tm\,/2 mw ./ m r, 1 „ .. . _ ^3 n .4 



-^[^(DJM^DJVifW^offdi agonal H ^ e4 < ( 176 ) 
m \ e e 4 ' 



1 v' Y u A u 

5 U lr = ^[LlV^(Q)A(U)V K (U)R u } A u =0 = " ^ | 

^. - ^MraW^mH 1 ^) e 4 

w m I g3 £ 6 £ 2 

r M(l)Wx(£)ii.] = ^| ^ ^ £ 2 

£ 5 £ 3 £ 2 

where the off diagonal parts are extracted except for stop mass matrix and 5^ L , 5 RR are omitted. 
4.4 Flavor and CP violation 

The off diagonal elements of sfermion mass matrices contribute to flavor and CP violation through the sfermion 
exchange, on which are imposed severe constraints. Based on the estimations of the flavor and CP viola- 
tions with the mass insertion approximation, the upper bounds for each elements are given in Table 4, where 
Mq = Af(gluino) = M(squark), Ml = Af(slepton) = A/(photino) are assumed [22]. Note that there is another 
suppression factor in 5^ R as 



V 
m 



(180) 



The most stringent bound for Ml is given by /i — > as 

/ M L 
V300GeV 



1< 1.5 x 10- 2 f-^- ) : M L > 2250GeV, (181) 



and the one for Mq is given by ex as 



( M Q 



e 2 5 = 3 x 10- 3 < 4.4 x 10~ 4 loooGeV J Mq > 6820GcV - ( 182 ) 

Note that if Q\,i were ^-singlets, then (<5^ £ )i2 would be O(l) and the most stringent bound for Mq would be 
given by 

1 < ^ X 10 " 3 ( lOb^ev) : ^> 156TcV - ( 183 ) 
Comparing Eq.(182) and Eq.(183), one can see that 54 softens the SUSY flavor problem very efficiently. 
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Before ending this section, we discuss the problem of a complex flavon VEV. If the relative phase of two VEVs 
(Di) , (D 2 ) exists, we must include 



K{D C 



[-£>?£> 2 +£>;£>i] (£>S)*(£>§) 
M 2 P 



+ h.c. 



in Kahlcr potential, then redefinition of supcrficlds arc modified as follows 



m 

D< 2 



= K'{D) 




1 e 3 
V K (D) = I e 3 1 



(184) 



(185) 



1 



Therefore the mass matrix and mixing matrix of down quark sector and off-diagonal matrix of squarks are 
modified as follows 



M' d 



R,i — 



r 

e 2 
1 



°RR — 



c 

r 2 





(186) 



As the result, the most stringent bound for Mq is changed into Mq > 68TcV. This suggests new mechanism is 
needed to suppress CP violation. We leave this problem for future work. 



upper bound 



observable 



parameter 



order 



coefficient 



0^)12(^)12 

(^L)l2 



,2.5 



5.6 x 1(T 3 
8.0 x 10~ 2 
8.8 x 10~ 3 



lOOOGeV 



(Amjc)exp 

3.49xlO- 12 AfeV 



Am B 



(^)l3(^)l3 

(*£i)l3 
(Ql3 



.1.5 



3.4 x 10~ 2 
1.9 x 10" 1 
6.3 x 10~ 2 



lOOOGeV 



(Am B )exp 

3.38xlQ- u, MeV 



Amu 



(<5L)i2(^ 

(^L)l2 
&12 



RR)i2 



1.1 X 10~ 2 

6.2 x 10~ 2 
1.9 x 10~ 2 



Mq 



lOOOGeV 



(Ami 



1.26xW- 11 MeV 



ImP£ L ) 12 (^) 12 ]| 



\MsE L )h\ 



,2.5 



llm^ ^ 2 



4.4 x 10~ 4 
6.4 x 10~ 3 
7.0 x 10- 4 



( m q 

^lOOOGeV 



2.24x10- 



LR1\1\ 



1.5 x 10~ 2 
3.4 x 10~ 6 



y ( Ml ' 
V300GeV/ 



(BR(^e~f)) e 



2.4x10- 



(€Z)l3 



r — > ej 



(8 E 



LR)13 



4.3 

1.6 x 10~ 2 



x ( Ml ' 

A V300GeV. 



(Bfl(r->e 7 )) e 
3.3xl0- s 



t — > fj/y 



(5 E 



L R .m 



1. 



4.9 

x 10~ 2 



/ Ml \ 2 / (BR(t^^)) cxp 
V 300GeV / V 4.4xl0- s 



3.1 x 10" B 
1.6 x 10~ 6 



M Q 



lOOOGeV 



)( 



(rfrQe 



2.9xl0~ 26 ecm 



M*£«)u| 



1.7 X 10~ 7 



V 300GeV - 



1.05xl0- 27 ecm 



Table 4: Experimental constraints for the off diagonal elements of sfermion mass matrices from meson mass 
splittings Amx,Ams,Amo, CP violating parameter e^, lepton flavor violations li — > Ijj and electric dipole 
moments of neutron d n and electron d e . The predictions of our model for each parameters are given in "order" 
column. The dependences of each upper bounds on experimental values are given in "coefficient" column. 



5 Cosmological Aspects 

Based on our model, we consider the scenario to reproduce the cosmological parameters given as follows [19] 

n ~ n A + n b + n C DM^i, (187) 
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n A = 0.73 ± 0.03. (188) 

Q b h 2 = 0.0225 ±0.0006, (189) 

VcDMh 2 = 0.112 ±0.006, (190) 

h = 0.704 ± 0.025. (191) 

For f2b, we adopt leptogenesis as the mechanism to generate baryon asymmetry. For Qcdm, we assume that 
dark matter consists of singlino dominated neutralino. 

5.1 Leptogenesis 

In general, leptogenesis scenario to generate baryon asymmetry causes over production of gravitino in supersym- 
metric model. This problem can be avoided in the case neutrino mass is generated by small VEV of neutrinophilic 
Higgs doublet [23]. 

In the diagonal RHN mass basis, superpotential of RHN is given by 

/ Y& Y» \ ( Nf\ 3 
W N = Y, e 3 Hf{L u L 2 ,L 3 )\ Y* yH 3 N% + - £ M a N c a N c a , (192) 

V Y» 2 Y» 3 J \N§ J 

where we assume accidental mass hierarchy as follows 

Mi = 10 3 - 5 , M 2 = M 3 = 10 4 (GcV). (193) 

Note that these particles are enough light to create in low reheating temperature such as 10 7 GeV without causing 
gravitino over production [9]. The interactions of right-handed sneutrinos (RHsNs) are given by 

o m 2 

-C N = Yl e 3 HV(L u L 2 ,L 3 ) | Y» 2 M 2 

i=l,2 

+ Y ^{hMM I y| 2 I I Ni |, (194) 

i=l,2 \ Y i 32 

where the contributions from A-terms are neglected. The Z 2 breaking scalar squared mass terms 

F* F 

K D B M | 2 " m)*NZ + ■■■] + h.c. = em 2 [(NcyNC + ■■■] + h.c (195) 

fill in the zeros of sneurino mass matrix and gives 

M 2 em 2 em 2 \ / e 

M 2 m 2 ~ Ml e 3 1 e 2 I , (196) 




M 3 > 




M 3 


1 ( ^ 


M 3 J 




Ni \ 








NI J 





em 2 m 2 M| ' 



where 0(e) suppressions of Z 2 breaking terms are assumed without any reason. Note that the 0(e 3 ) elements 
are originated from small Z 2 breaking parameters and small Y® as 



2 

<S>\2 3 



e(Y*y ~ e 3 . 

In the diagonal RHsN mass basis, the interaction terms given in Eq.(194) are rewritten by 

iVf 

~£n = ^ f 3 H^ r (Li,L 2 ,L 3 

i=l,2 

Y,^{hM,h)[ ^li Y» 2 Y» 23 N° ]. (198) 

i=l,2 
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The lightest RHN n\ does not receive above corrections and remains decoupled. Therefore lepton asymmetry is 
generated by the out of equilibrium decay of the lightest RHsN iVf . 

Following [24], the CP asymmetry of sneutrino iV-j 3 decay is calculated as follows 



ei 



i;E4^^)' ( 199 ) 



4tt ^ ifn 



= ^inl+^ + M-, (200) 
x x — 1 

M 2 

^ = ^ (20i) 



K » = E E^MiJ^Wi)*- ( 202 ) 

h=l,2 1=1 

From the naive power counting, we obtain 

A'n~e 12 , K 12 ~ A' 13 ~ e 9 , e x ~ e 6 . (203) 
Using ei , the B — L asymmetry generated via thermal leptogenesis is expressed as 

— {B — L) f = k—, g* =341.25, (204) 
9* 

where g* is the total number of relativistic degrees of freedom contributing to the energy density of the universe 
and dilution factor k is defined as follows 

« ~ (205, 

K = W (206) 
2H(M 1 ) y ' 

r(M0 = (207) 

87T 



"< M >> - <208) 

By the EW sphalcron processes, the B — L asymmetry is transferred to a B asymmetry as 

24 + ANh 1 

where Nh is number of Higgs doublets which are in equilibrium through Yukawa interactions, for example Nh = 1 
for SM and Nh = 2 for MSSM. In any way A^-dependence is not important for our rough estimation. For our 
parameter values, we obtain K ~ O(l) and 

%-10~ 10 , (210) 

which is consistent with observed value 

r] B = 7.04% = 6.1 x 10~ 10 . (211) 

Requiring the effective interaction 

is decoupled in order to avoid too strong wash out, we impose the condition as follows 



e 12 T 3 7r 2 a*T 4 
r~ , _ <H = J 3 *\ , (213) 

87T 3 M 2 2 V 90M p 
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which gives upper bound for temperature as 

T < 10 4 GeV. (214) 

This condition is always satisfied after the decay of JVf starts. 

This scenario is different from conventional one in the point that neutrino mass 

^~^~(^)\leV~O(0.01eV), (215) 

is realized by the small VEV v u — O(lGeV). 

5.2 Dark matter 

Here we calculate the relic abundance of LSP which corresponds to singlino dominated neutralino in our model 
[25]. The most dominant contribution to annihilation cross section of LSP is given by the interaction with Z 
boson. If the mass matrix given in Eq.(97) is diagonalized by the field redefinition as 

(hf)° \ ( V a * * \ / X°i \ 

(/if) = H * * Xi, 2 , m x o <m x o 2 <m x o 3 , (i = 1,2), (216) 

*i J \V C * * J \ Xi,3 J 
the interaction with Z boson is given by 

= Gixi^Z^xl.+iGU^frz^PLf + tGlMfrz.Pnf, (217) 



where 



G(xh) = \ {Va V?)y/fi+gl = 0.372(|K| 2 - \V b \ 2 ) (218) 
G(e L ) = 0.200, G{e R ) = -0.172, G{v L ) = -0.372, 

G{u L ) = -0.257, G(u R ) = 0.115, G(d L ) = 0.314, G{d R ) = -0.057, (219) 



Oi Y (m z ) = 0.0101687, a 2 {m z ) = 0.0338098 (220) 



are used. 

The formula for the relic abundance of cold dark matter is given by 

n CDM h> = 8J6 x m- 11 *^, (221 ) 

(a + 3b/x F )GeV 2 1 ; 

0.0955mpm v o(a + 66/ xf) 
x F = In '—, (222) 

{g*XF) 2 

rap = 1.22 x 10 19 GeV, 

.g* = 72.25 (T F = m x o/x F < m T ), (223) 



1 2 



4m^ — m z 



(224) 



6 = 



£ s c2 <*> (s^r^i) 2 [c2(/t) + G2(ft)1 - r <225) 

where rrif <C ti x o is assumed. Substituting the values given in Eq.(218) and Eq.(219) and following values 

2m y o 

m z = 91.1876, m 6 = 4.18, m T = 1.777 (GcV), r= ^ (226) 



in Eq.(224) and Eq.(225), we get 



2 



,- 8 | G(xji) 
y r 2 - 1 J 

2 / . „ . \ 2 



1.77 x 10~ 8 —fzT ( GeV )' ( 22? ) 



6.436 x 10- 6 j — 0.013 x 10~ 6 I (GcV^), (228) 
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from which the formula is rewritten as follows 



%F 



^CDMh 2 



= In 



0.01 



0.0177 + — (6.436r 2 - 0.013) 
0.10306a; f 



6.25 x 10 s 

(r 2 - l) 2 



2 H^f), 



0.01 



0.0177 + ^(6.436r 2 - 0.013) 1 -r^ 



(229) 
(230) 



where quark and lepton masses are neglected except for bottom and r. Since the two LSPs Xi i > X2 1 have the 
same mass and the same interactions, they have the same relic abundance. Therefore the required relic abundance 
of one LSP is QcDAth 2 = 0.055. For the allowed range given in Eq.(104), the required values for A4 ; 5 to reproduce 
observed relic abundance of dark matter are given in Table 5. The allowed ranges for A 4i 5 are very small. Note 
that we should not impose LEP bound (m x o > 46GeV) on this LSP, because Z — s- Xi,iXi,i 1S strongly suppressed 
by the factor (|Vq| 2 — |Vfc| 2 ) 2 /2 ~ 0.005 and the contribution to invisible decay width is negligible as follows 



T(Z -> xliX°i,i) ~ (0.6 x 2/3)0.005r(invisible) ~ l.OMeV, 
r(invisible) = 499.0 ± 1.5MeV[19], 



(231) 
(232) 



where phase space suppression factor ~ 0.6 and the ratio of LSP number and neutrino number 2/3 are multiplied. 



A 4 


A 5 


m ± 


*i,i 


TO v 

-M,2 




V a 


v b 


Xf 


T F 




0.44 


0.57 


141.2 


36.47 


142.38 


178.84 


0.3667 


0.2225 


22.90 


1.592 


0.0552 


0.42 


0.56 


130.1 


36.52 


131.09 


167.61 


0.3721 


0.2321 


22.90 


1.595 


0.0550 


0.40 


0.55 


119.1 


36.61 


119.91 


156.52 


0.3776 


0.2429 


22.90 


1.599 


0.0551 


0.38 


0.54 


108.0 


36.79 


108.63 


145.42 


0.3836 


0.2554 


22.91 


1.606 


0.0549 


0.37 


0.53 


102.5 


36.60 


103.10 


139.69 


0.3885 


0.2591 


22.91 


1.597 


0.0550 



Table 5: The parameter sets (A4, A5, mi = ^i v ' s ) which reproduce observed relic abundance of dark matter. The 
dimensionful values are expressed in GeV units. 



5.3 Constraint for long-lived massive particles 

Finally we consider long-lived massive particles which are included in our model, G-Higgs, Savons and the lightest 
RHN. Such particles are imposed on strong constraints from cosmological observations. 
The superpotential of G-Higgs sector gives degenerated G-higgsino mass as 



M g = fc<diag(l, 1, 1), 
which receives S4 breaking perturbation from Kahler potential given by 



K(G) 



\G a \ 



Ml i 



2D 2 G 1 \ 2 



-VaDi - D 2 )G 2 \ 2 + |(VaDi - D 2 )G 3 \ 2 + (G -> G c ) 



= \G a \ 2 + J2^e 2 \G a \ 2 + (G^G c ), 



(233) 



(234) 



which solves the mass degeneracy, however generation mixing is not induced. Neglecting 0(e 2 ) corrections and 
contributions from D-tcrms except for the contribution from S3, the G-Higgs mass terms are given by 



-C D m 2 \G a 



\G c a \ 2 - [kA k S 3 G a G c a + h.c] + \kS 3 G a \ 2 + \kS 3 G c a 



tUttD\ 



+ \kG a G c a + X 3 H^H 3 
from which we obtain three same 2x2 matrices as 
M 2 (G) = 



\g 2 x [b\S 3 \ 2 -2\G a \ 2 -Z\Gl\ 2 } 2 . 



m 2 G + {kvi s ) 2 ~l0g 2 M) 2 



X 3 kv' u v' 



hkv' u v' d - kA k v' s 



(kv>) 2 - 15g 2 x (v' s ) 2 



(235) 



(236) 
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The mass spectrum of G-Higgs and G-higgsino is given in Table 7 and the lightest particle of them is lighter 
G-Higgs scalar G_ . The dominant contributions to the G_ decay are given by the superpotential 



M 2-— — » Ml 

-^=Y QQ Q 3 Q 3 (Gi + G 2 + G 3 ) + -j=Y UE U§Ei(Gi + G 2 + G 3 ), (237) 
Y QQ = y ue = (<M\ \ 2x 10 _u (238) 



from which we obtain 



C G = ^A^Y UE (e c 3 u c 3 + u c 3 e c 3 )G 1 + ^=A^Y^(2u 3 d 3 + 2d 3 u 3 )G 1 . (239) 



For simplicity, we assume G ~ G_ then decay width of G_ is given by 

M(G_) 



r(G_) 



(Y QQ ) 2 , (240) 



16tt 

where the rcnormalization factors are calculated based on the RGEs given in appendix A as follows 



Al E F = J auE[Ms \ =4.9, A Q S = J ggg^4 = i2.8. (241) 
* F \J a UE (M P ) RF \ja QQ {M P ) K ' 

Substituting these values in Eq.(240) we obtain the life time of G_ as 

r(G_) = f ^- y= 3.8xl0- 29 (^iy 1 (r^)- 2 sec. (242) 

Since the existence of a particle which has longer life time than 0.1 second spoils the success of BBN [9], we must 
require r(G_) < O.lsec which impose constraint as 



1.9 x 10" 14 



2 



M(G-) > ^ • yQQ ) TeV - (243) 

The G-Higgs exchange may contribute to proton decay, however it seems that the suppression of power of e is 
too strong to observe proton decay [26]. 

The five of six flavon multiplets ^ a ,^ have O(lTeV) masses which are enough small to product them 
non-thcrmally through the U(l)z gauge interaction. The lightest flavon (LF) is quasi-stable and should not 
produced so much in order not to dominate SIcdm ■ Solving the Boltzmann equation with the boundary condition 
tilf(Tjih) — 0, we get relic abundance of LF as 1 

tl LF h 2 = 2.0 x 10- 8 ( -^-V ( H^XV = 2.0 x 10- 6 ( -^-V [27]. (244) 

Requiring the LF does not dominate dark matter as Q^ph 2 < 0.01, the upper bound for reheating temperature 
is given by 

T R < 10 6 GeV, (245) 

which is consistent with our leptogenesis scenario. 

The life time of LF is estimated as follows. The LF can decay, for example through the operator 



W = — ^ ^ l — U- 1 - 2- , (246) 

2(V + <S>) 2 2V 2 v v; 



^ince the U(l)z charge of <E> in [27] is two times larger than one in the present model, we multiply the equation for fl^ph 2 given 
in [27] by the factor 2 2 . 
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the decay width and life time are given by 



T(LF -> IIHH) = ^ ( ^}f}' I O(O.l) ~ 10- 29 eV, 



16tt \32ir 2 M 2 V t 

t(LF) ~ 10 14 scc- 10 7 years[27], (248) 

which suggests LF does not exist in present universe. Note that three and two body decays are suppressed by 
small VEV v u . 

The lightest RHN n\ behaves like LF because there is no distinction between N c and <I> C under the gauge 
symmetry. Integrating out N£ and Xz in the Lagrangian 

C D gz {niXz(N^)* + ipXz**) + e e N?lhY, (249) 

where ("0, V J') means super-multiplet and some factors are omitted for simplicity, we get 

C eS=J^^(nimhY)% (250) 
from which the life time of n\ is given by 

Y[n\ -> ip^lh u ) = ™l ( g ^2 2 ") ~ ( ^-) 5 T(LF ->■ IIHH) ~ lQ- 22 eV, (251) 



167r(327r 2 ) 2 \{gzV)Ml) \M S 

Ms 
Mi 



r K) ~ (tt t ( LF ) ~ 10 12 sec - 10 5 years. (252) 



6 Conclusion 

In this paper we consider S4 flavor symmetric extra U(l) model and obtain following results. 

• With the assignment of flavor representation to reproduce quark and lepton mass hierarchies and mixing 
matrices, SUSY flavor problem is softened. 

• Proton decay through G-Higgs exchange is suppressed by flavor symmetry. 

• Observed Higgs mass 125 — 126GeV is realized with stop lighter than 2TeV which is within the testable 
range in LHC at y/s = 14TeV. 

• The partial gauge coupling unification at Mp is realized by adding 4-th generation Higgs and left-handed 
lepton which play the role to break U(l)z gauge symmetry. 

• The allowed region for lightest chargino mass is given by 100 — 140GcV when we assume LSP is lightest 
singlino dominated ncutralino. 

• The extra Higgs doublets play the role of neutrinophilic Higgs which is needed for low temperature lepto- 
gengesis without causing gravitino over production. 

• The shorter life time than 0.1 second of G-Higgs is realized. 

• The over productions of flavon and lightest RHN are also avoided. 
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A RGEs 



O(l) coupling constants of our model consist of gauge coupling constants and trilinear coupling constants defined 
by 



W 



D \3S3H3 + \±H 3 {S\Hi + S2H2) + XsiSiHi 

rU TjU, 



kS^dGl + G 2 G% + G3G3) + Y 3 U H£Q 3 U£, 



(253) 



from which the fine structure constants are defined as follows 



ay 



_9± 
4tt' 

47T 



«2 



4tt' 



1 :S 



Oth 



\2 

4tt' 



4tt' 
Q4 = 



ax 



\2 

A 4 



4tt 



05 



4tt' 
= A| 

47T 



OA' 



4tt' 
47r' 



We define the step functions as follows 
9(x)-- 



vi = 
Mi 



1 x > 
x < 

-- 10 115 GeV, 



4 = 0(M-A/L 4 ), #5 = 

M L4 = 2.2 x 10 14 GcV, 



0(p-Ml 5 ), 
M L . = 2.4 x 10 17 GcV. 
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The beta functions are given by 
. day 
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50 
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80 03 
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where we include only the contributions from 012,3 in 2-loop order terms. The RGEs for trilinear coupling 
constants are given by 



at 6a 



16 o 13 1 5 n 

an + 2a4 — a 3 - 3a 2 — —ay - -ax - -^a z 0i 

3 9 2 o 

19 5 

oth [ 3a t + 4a^ + 2a 4 + 2a 5 + 9afc - 3a 2 - ay — —a x - -az@i 

6 6 



= «4 3at + ah + 5a4 + 2a5 — 3a 2 



(27T) 
(27T) 



das 

yr 



C*5 I OLh + 2tt4 + 5a 5 — 3a 2 — Oty 



ay 
19 



19 5 
6 6 



aft 2a/i + Haft 



16 



-a 3 



-ay 



6 
19 



a x 



3 9 6 
We define gaugino mass parameters and A-parameters as follows 

C D 



a x 



a z Vi 



a z 6, 
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(267) 



-My Ay Ay - iA/ 2 A 2 A 2 - ^M 3 Af Af - -MxXx^x - \m z \ z \ z 



\ 3 A 3 S 3 H%Hg + XiAiH^S^ + S2H2) + XsMSiHi + SzH$)Hg 



+ kAkS^GxGl + G 2 G L 2 + G3G3) 



Y 3 u A t H^Q 3 U^ + h.c. 



(268) 
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The RGEs for gaugino mass parameters are given by 

dM Y 



(27T). 
(27T) 
(27T) 
(27T) 
(27T) 



dt 

dM 2 
dt 
dM 3 

dt 
dM x 
dt 
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= ay My 
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= a z M z 
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65 20 

1 
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(269) 
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where we take account of 2-loop contributions only for M 3 . The RGEs for A-parameters are given by 
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RGEs for scalar squared masses are given by 
dm%. 
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where 



(2^)—^ - a fc M fc 2 - -a 3 M 2 - -a r M Y - -a x M 2 - ^a z M 2 e T , (288) 

Ml = Aj + m% 3 + m 2 H u + m 2 , 
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= A 2 
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= m G 2 : 


= m% 3 = 



Ml =A 2 + mj; i + ^ + ; 



(289) 



2 _ 2 _ 2 2 
777' it u — Tflirn — 777 it n. TYlri 



(290) 
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are held. At /i = M/, we add U(l)z D-term corrections as follows [28] 

m 2 x {Mi-0) = m 2 x {Mi + 0) + Am 2 x (291) 
A?77q = Amfjc = Am%c = Am 2 HD = Am G c = Airig = — tt^dt^ 

Am 2 D c = Am\ = Am 2 HU = Am G = --m 2 DT , (292) 

m 2 DT = ITeV 2 > 0. (293) 

We solve these RGEs using following boundary conditions. At SUSY breaking scale (fi = Ms = ITeV), we put 
by hand as follows 

A 3 = 0.37, A 4 = 0.4, A 5 = 0.55, Y t = Y 3 U = 1.0, k = 0.5, M 3 = lOOOGeV, M Y = 200GcV, 
m 2 Qi = 3.00, mljc = 1.00, m 2 HU = m 2 HD = m 2 s = 2.00, m 2 G = 5.50, m 2 GC = 7.00 (TeV 2 ), 

m 2 H u H o S3 -> Eq.(67)(68)(69). (294) 
At reduced Planck scale (/j = Mp = 2.4 x 10 18 GeV),we put by hand as follows 

a 2 =a 3 = 0.125, a x = a z = a Y = 0.209, M 2 = M 3 , M x = My = M z , A t , 3 ,4,5,k = 0, 

m 2 La = m 2 EC = m 2 D c = mfjc = m,Q = 0. (295) 

Note that gauge coupling constants do not satisfy the conventional unification as 

3 . . 

ay = -a 2 3 . (296) 
5 

The renormalization factors of first and second generation Yukawa coupling constants Y u ' d ' e and single G-Higgs 
coupling constants defined by 

W G = Y QQ Q 3 Q 3 (G\ + G 2 + G 3 ) + Y UE U^E C 3 (G 1 + G 2 + G 3 ), (297) 

are given by 



\ r A ?f?\ > «A = ^, A = u,d,e,QQ,UE, (298) 
y a A {M P ) 47r 

which are calculated by RGEs as follows 

1 da u 16 13 1 5 

(2tt) — = 3a t + a h + 2a4- —a 3 -3a 2 - —a Y ~ -a x - -a z 0!, (299) 

a u at 3 9 2 6 

1 dot 16 7 7 5 

(2ir) — = a h + 2a 5 - —a 3 - 3a 2 - -ay - -a x - -a z 6i, (300) 

ad at 3 9 6 6 

(2tt) — = ah + 2a<5 - 3a 2 - 3cty - \a x - ^a z 6i, (301) 
a e at 6 6 

(0 \ 1 dauE O a. 16 28 1 5 ft (W>\ 

(2vr) -r— = 2a t + a k ~ —a 3 - -ay - -a x - -a z tii, (302) 

ajjE at 3 9 2 6 

(27t)—!— da QQ _ 2a t + a k ~ 8a 3 - 3a 2 - \a Y - \a x - %a z R i, (303) 
a QQ at 3 2 6 

where the contributions from $ a , $3, Di are neglected. The results are given in Table 6. 
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parameter 


fi = M s (Mj) 


/x = M P 


parameter 


H = M s 


[i = M P 


ay 


0.010442 


0.209 




3.00 


0.9912 


a 2 


0.032482 


0.125 




1.00 


3.1451 




0.089430 


0.125 




-0.1723 


16.4770 


ax 


0.010552 


0.209 


P 

H u 


2.00 


2.4671 


a z 


(0.015162) 
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m' 1 D 
£ 


2.6811 


6.2912 


ctt 


0.079577 


0.006455 




2.00 


1.6459 


ah 


0.010894 


0.016086 


03 


-2.2105 


10.4216 


Qf4 


0.012732 


0.011761 


2 

m S 


2.00 


6.4914 


"5 


0.024072 


0.011309 


m G 


5.50 


1.6852 


Otk 


0.019894 


0.001014 


m z Gc 


7.00 


2.2501 


M Y 


0.2 


3.68582 


m Q 


5.9677 


0.0 


M 2 


0.49889 


1.68366 


"2* 

mfja 


5.7726 


0.0 


M 3 


1.0 


1.68366 


nine 


4.8538 


0.0 


M x 


0.20228 


3.68582 


mi 


1.0559 


0.0 


M z 


(0.28352) 


3.68582 


r 

n 


1.7796 


0.0 


A t 


-1.78127 


0.0 


a u 


0.259802 


0.01 


A 3 


1.51978 


0.0 


ad 


0.523940 


0.01 


A, 


0.40541 


0.0 


a e 


0.037691 


0.01 


A, 


-0.93253 


0.0 


auE 


0.238230 


0.01 


A k 


-3.05177 


0.0 


a QQ 


1.638921 


0.01 



Table 6: Each boundary values of the solutions of RGEs. The dimensionful parameters are expressed in TeV units. 
The experimental values of gauge coupling constants give ay (Ms) = 0.010445, a 2 (Ms) = 0.032484, a3(Ms) = 
0.089514 which are calculated based on SM RGEs. The values of az,Mz at low energy side are given by the 
values at fi = Mi in brackets. 



B The multiplication rules of S4 

The representations of S4 are 1, 1', 2, 3, 3' [29]. Their products are expanded as follows. 

' X2V3 ~ 2=32/2 

2=32/1 - 2=12/3 
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£2 2/3 + 3:32/2 
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xiy 2 + 2=22/1 
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2:22/3 - 2=32/2 \ / x 2 y 3 + x 3 y 2 

2=32/1-2=12/3 + 2; 3 yi+xiy 3 | , (305) 
2=i2/2-a=22/i / 3 \ xiy 2 +x 2 yi , 3 , 

2x22/1 \ / 2xiyi 

-v / 3xij/2 - x 2 y 2 + -X12/2 + \/3x 2 y 2 | (306) 

V3xiy 3 - x 2 y 3 / 3(3/) \ -xiy 3 - V3x 2 y 3 / 3 
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x (vh 



0)i' x {y)v 



-X2V 

xiy 



(309) 
(310) 



C Mass bounds of new particles 



particle 


mass 


exp 


particle 


mass 


exp 


if (lightest even) 


125.7 


125 - 126 [2] 


xf 


1486 


> 295 [16] 


T+ 


1882 


> 560[16] 




493 


> 295 [16] 


T_ 


1178 


> 560[16] 


xS 


199 


> 46 [19] 


G+ 


3908 




X°2 


493 


> 62.4[19] 


G- 


1737 


(> 683) [16] 


X b 3 


1481 


> 99.9[19] 


Ql,2 


> 2532 


> 1380[16] 


xi 


1487 


> 116 [19] 




> 2493 


> 1380[16] 


xt 


2004 




u l,2,3 


> 2395 


> 1380[16] 


xi 


2208 




Ll,2,3 


> 1393 


> 195[16] 


xt 


119 


100- 140 [15] [16] 


^1,2,3 


> 1490 


> 195[16] 


xl 


36.6 




Hi 2 {even, odd, ±) 


1056 


> 93.4[19] 


xh 


120 


> 116 [19] 


Hi 2 i even i odd, ±) 


821 


> 93.4[19] 


xl 


157 


> 116 [19] 


Si^even, odd) 


2052 




9 


2000 




H^(even, odd, ±) 


2279 


> 93.4[19] 


\ 9 

A 3 


1000 


> 1000[16] 


S3(even) 


2102 




Z' 


2102 


> 1520 [7] 



Table 7: Mass values of new particles calculated based on our assumption and corresponding experimental 
constraints in GcV units. The capital letters means bosons and the Greek characters and the small letter 
mean fermions. The equations which are used to calculate mass values, are Eq. (10), (71), (74), (75), (76), (77), 
(83), (85), (95), (96), (97), (98). Each equalities in "mass" column correspond to imposing the boundary conditions 
as m 2 x = 0(X = Q,Uf,D c a ,L a ,E^) at /x = Mp. We adopt the mass bound for stable stop as one for lighter 
G-Higgs (G-) in bracket, under the assumption that G_ is lighter than g and G+. We adopt the mass bound 
for CP-odd Higgs boson in supersymmetric model as ones for extra Higgs bosons (H^ 2 i H3). 

The mass bound of the lightest chargino (xf ) is given by 3-lepton emission through EW direct process 
X1X2 ~^ W^ZxiXi- This neutralino \i corresponds to Xii m our model. Under the assumption that slepton 
decouples and LSP (x?) is massless, excluded region of chargino mass is given by 140 < M(xf) < 295GeV [16], 
or M(xf) < 330GeV [17]. These constraints are not imposed on the chargino in the mass range 

M( X f) = M(x° 2 ) < Mix !) + m z . (311) 

In this case Z and following two lepton emissions are suppressed. Taking account of LEP bound M(xf) > lOOGeV 
[15], we consider the allowed region given by 

100 < M(xf) < 140 (GeV). (312) 
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